Trabajo 1 TAE - Salud
Librerías
library(tidyverse)
library(data.table)
library(readxl)
library(GGally)
library(leaflet)
library(raster)
library(mltools)Lectura de datos
raw_encuesta <- fread("databases/calidad_vida_ok.csv", encoding = "UTF-8") %>%
as_tibble()
raw_encuesta$encuesta_calidad.barrio %>% unique() %>% sort()## [1] "AGUAS FRÍAS"
## [2] "ALDEA PABLO VI"
## [3] "ALEJANDRÍA"
## [4] "ALEJANDRO ECHAVARRÍA"
## [5] "ALFONSO LÓPEZ"
## [6] "ALTAMIRA"
## [7] "ALTAVISTA"
## [8] "ALTAVISTA CENTRAL"
## [9] "ALTOS DEL POBLADO"
## [10] "ANDALUCIA"
## [11] "ANTONIO NARIÑO"
## [12] "ARANJUEZ"
## [13] "AREA DE EXPANCION SAN CRISTOBAL"
## [14] "ÁREA DE EXPANSIÓN ALTAVISTA"
## [15] "ÁREA DE EXPANSIÓN BELÉN RINCÓN"
## [16] "AREA EXPANSION PAJARITO"
## [17] "AREA EXPANSION SAN ANTONIO DE PRADO"
## [18] "ASTORGA"
## [19] "AURES Nº 1"
## [20] "AURES Nº 2"
## [21] "BARRIO CAYCEDO"
## [22] "BARRIO COLOMBIA"
## [23] "BARRIO COLON"
## [24] "BARRIO CRISTÓBAL"
## [25] "BARRIO FACULTAD DE MINAS"
## [26] "BARRIOS DE JESÚS"
## [27] "BARRO BLANCO"
## [28] "BASILIA"
## [29] "BATALLÓN GIRARDOT"
## [30] "BELALCAZAR"
## [31] "BELÉN"
## [32] "BELENCITO"
## [33] "BELLO HORIZONTE"
## [34] "BERLÍN"
## [35] "BERMEJAL- LOS ÁLAMOS"
## [36] "BETANIA"
## [37] "BLANQUIZAL"
## [38] "BOLIVARIANA"
## [39] "BOMBONA Nº 1"
## [40] "BOMBONA Nº 2"
## [41] "BOQUERÓN"
## [42] "BOSQUES DE SAN PABLO"
## [43] "BOSTON"
## [44] "BOYACÁ"
## [45] "BUENOS AIRES"
## [46] "BUGA PATIO BONITO"
## [47] "CABECERA ALTAVISTA"
## [48] "CABECERA SAN ANT DE PR."
## [49] "CABECERA SAN CRISTÓBAL"
## [50] "CALAZANS"
## [51] "CALAZANS PARTE ALTA"
## [52] "CAMPO ALEGRE"
## [53] "CAMPO AMOR"
## [54] "CAMPO VALDÉS Nº 1"
## [55] "CAMPO VALDES Nº 2"
## [56] "CARIBE"
## [57] "CARLOS E RESTREPO"
## [58] "CARPINELO"
## [59] "CASTILLA"
## [60] "CASTROPOL"
## [61] "CATALUÑA"
## [62] "CERRO EL VOLADOR"
## [63] "CERRO NUTIBARA"
## [64] "CIUDADELA NUEVO OCCIDENTE"
## [65] "CORAZÓN DE JESÚS"
## [66] "CÓRDOBA"
## [67] "CORREGIMIENTO PALMITAS"
## [68] "CRISTO REY"
## [69] "CUARTA BRIGADA"
## [70] "CUCARACHO"
## [71] "DESCONOCIDO"
## [72] "DIEGO ECHAVARRÍA"
## [73] "DOCE DE OCTUBRE Nº 1"
## [74] "DOCE DE OCTUBRE Nº 2"
## [75] "EDUARDO SANTOS"
## [76] "EL ASTILLERO"
## [77] "EL CARMELO"
## [78] "EL CASTILLO"
## [79] "EL CERRO"
## [80] "EL CHAGUALO"
## [81] "EL COMPROMISO"
## [82] "EL CORAZÓN"
## [83] "EL CORAZÓN EL MORRO"
## [84] "EL DANUBIO"
## [85] "EL DIAMANTE"
## [86] "EL DIAMANTE Nº 2"
## [87] "EL JARDÍN"
## [88] "EL LLANO"
## [89] "EL NOGAL - LOS ALMENDROS"
## [90] "EL PATIO"
## [91] "EL PESEBRE"
## [92] "EL PICACHO"
## [93] "EL PINAL"
## [94] "EL PLACER"
## [95] "EL PLAN"
## [96] "EL PLAYON DE LOS COMUNEROS"
## [97] "EL POBLADO"
## [98] "EL POMAR"
## [99] "EL RAIZAL"
## [100] "EL RINCÓN"
## [101] "EL SALADO"
## [102] "EL SALVADOR"
## [103] "EL SOCORRO / LA GABRIELA"
## [104] "EL TESORO"
## [105] "EL TRIUNFO"
## [106] "EL UVITO"
## [107] "EL VELÓDROMO"
## [108] "EL YOLOMBO"
## [109] "ENCISO"
## [110] "ESTACIÓN VILLA"
## [111] "ESTADIO"
## [112] "FÁTIMA"
## [113] "FERRINI"
## [114] "FLORENCIA"
## [115] "FLORIDA NUEVA"
## [116] "FRANCISCO ANTONIO ZEA"
## [117] "FUENTE CLARA"
## [118] "GERONA"
## [119] "GIRARDOT"
## [120] "GRANADA"
## [121] "GRANIZAL"
## [122] "GUAYABAL"
## [123] "GUAYAQUIL"
## [124] "HÉCTOR ABAD GÓMEZ"
## [125] "JESÚS NAZARENO"
## [126] "JUAN PABLO II"
## [127] "JUAN XXIII - LA QUIEBRA"
## [128] "KENNEDY"
## [129] "LA AGUACATALA"
## [130] "LA ALDEA"
## [131] "LA AMÉRICA"
## [132] "LA ASOMADERA Nº 1"
## [133] "LA ASOMADERA Nº 2"
## [134] "LA ASOMADERA Nº 3"
## [135] "LA AVANZADA"
## [136] "LA CANDELARIA"
## [137] "LA CASTELLANA"
## [138] "LA COLINA"
## [139] "LA CRUZ"
## [140] "LA CUCHILLA"
## [141] "LA ESPERANZA"
## [142] "LA ESPERANZA Nº 2"
## [143] "LA FLORESTA"
## [144] "LA FLORIDA"
## [145] "LA FRANCIA"
## [146] "LA FRISOLA"
## [147] "LA FRONTERA"
## [148] "LA GLORIA"
## [149] "LA HONDONADA"
## [150] "LA ILUSIÓN"
## [151] "LA ISLA"
## [152] "LA LADERA"
## [153] "LA LIBERTAD"
## [154] "LA LOMA"
## [155] "LA LOMA DE LOS BERNAL"
## [156] "LA MANSIÓN"
## [157] "LA MILAGROSA"
## [158] "LA MOTA"
## [159] "LA PALMA"
## [160] "LA PILARICA"
## [161] "LA PIÑUELA"
## [162] "LA PRADERA"
## [163] "LA ROSA"
## [164] "LA SALLE"
## [165] "LA SIERRA"
## [166] "LA SUCIA"
## [167] "LA SUIZA"
## [168] "LA VERDE"
## [169] "LALINDE"
## [170] "LAS ACACIAS"
## [171] "LAS BRISAS"
## [172] "LAS ESMERALDAS"
## [173] "LAS ESTANCIAS"
## [174] "LAS GRANJAS"
## [175] "LAS INDEPENDENCIAS"
## [176] "LAS LOMAS Nº 1"
## [177] "LAS LOMAS Nº 2"
## [178] "LAS MERCEDES"
## [179] "LAS PALMAS"
## [180] "LAS PLAYAS"
## [181] "LAS VIOLETAS"
## [182] "LAURELES"
## [183] "LLANADITAS"
## [184] "LÓPEZ DE MESA"
## [185] "LORENA"
## [186] "LORETO"
## [187] "LOS ALCÁZARES"
## [188] "LOS ALPES"
## [189] "LOS ÁNGELES"
## [190] "LOS BALSOS Nº 1"
## [191] "LOS BALSOS Nº 2"
## [192] "LOS CERROS EL VERJEL"
## [193] "LOS COLORES"
## [194] "LOS CONQUISTADORES"
## [195] "LOS MANGOS"
## [196] "LOS NARANJOS"
## [197] "LOS PINOS"
## [198] "MANILA"
## [199] "MANRIQUE CENTRAL Nº 1"
## [200] "MANRIQUE CENTRAL Nº 2"
## [201] "MANRIQUE ORIENTAL"
## [202] "MARIA CANO - CARAMBOLAS"
## [203] "MAZO"
## [204] "MEDIA LUNA"
## [205] "METROPOLITANO"
## [206] "MIRA FLORES"
## [207] "MIRADOR DEL DOCE"
## [208] "MIRANDA"
## [209] "MIRAVALLE"
## [210] "MONTAÑITA"
## [211] "MONTECLARO"
## [212] "MORAVIA"
## [213] "MOSCU Nº 1"
## [214] "MOSCU Nº 2"
## [215] "NARANJAL"
## [216] "NUEVA VILLA DE ABURRA"
## [217] "NUEVA VILLA DE LA IGUANA"
## [218] "NUEVOS CONQUISTADORES"
## [219] "OCHO DE MARZO"
## [220] "OLAYA HERRERA"
## [221] "ORIENTE"
## [222] "PABLO VI"
## [223] "PAJARITO"
## [224] "PALENQUE"
## [225] "PALERMO"
## [226] "PATIO BONITO"
## [227] "PEDREGAL"
## [228] "PEDREGAL ALTO"
## [229] "PICACHITO"
## [230] "PICACHO"
## [231] "PIEDRA GORDA"
## [232] "PIEDRAS BLANCAS"
## [233] "POPULAR"
## [234] "POTRERA MISERENGA"
## [235] "POTRERITO"
## [236] "PRADO"
## [237] "PROGRESO"
## [238] "ROBLEDO"
## [239] "ROSALES"
## [240] "SAN ANTONIO"
## [241] "SAN BENITO"
## [242] "SAN BERNARDO"
## [243] "SAN DIEGO"
## [244] "SAN GERMÁN"
## [245] "SAN ISIDRO"
## [246] "SAN JAVIER Nº 1"
## [247] "SAN JAVIER Nº 2"
## [248] "SAN JOAQUÍN"
## [249] "SAN JOSÉ"
## [250] "SAN JOSE DE LA MONTAÑA"
## [251] "SAN JOSÉ DEL MANZANILLO"
## [252] "SAN JOSE LA CIMA Nº 2"
## [253] "SAN JOSELA CIMA Nº 1"
## [254] "SAN LUCAS"
## [255] "SAN MARTÍN DE PORRES"
## [256] "SAN MIGUEL"
## [257] "SAN PABLO"
## [258] "SAN PEDRO"
## [259] "SANTA CRUZ"
## [260] "SANTA ELENA SECTOR CENTRAL"
## [261] "SANTA FE"
## [262] "SANTA INES"
## [263] "SANTA LUCIA"
## [264] "SANTA MARGARITA"
## [265] "SANTA MARIA DE LOS ÁNGELES"
## [266] "SANTA MÓNICA"
## [267] "SANTA ROSA DE LIMA"
## [268] "SANTA TERESITA"
## [269] "SANTANDER"
## [270] "SANTO DOMINGO SABIO Nº 1"
## [271] "SANTO DOMINGO SABIO Nº 2"
## [272] "SECTOR CENTRAL"
## [273] "SEVILLA"
## [274] "SIMÓN BOLIVAR"
## [275] "SUCRE"
## [276] "SURAMERICANA"
## [277] "TEJELO"
## [278] "TOSCANA"
## [279] "TRAVESÍAS"
## [280] "TRECE DE NOVIEMBRE"
## [281] "TRICENTENARIO"
## [282] "TRINIDAD"
## [283] "URQUITA"
## [284] "VEINTE DE JULIO"
## [285] "VERSALLES Nº 1"
## [286] "VERSALLES Nº 2"
## [287] "VILLA CARLOTA"
## [288] "VILLA DEL SOCORRO"
## [289] "VILLA FLORA"
## [290] "VILLA GUADALUPE"
## [291] "VILLA HERMOSA"
## [292] "VILLA LILIAM"
## [293] "VILLA NIZA"
## [294] "VILLA NUEVA"
## [295] "VILLA TINA"
## [296] "VILLA TURBAY"
## [297] "VOLCANA GUAYABAL"
## [298] "YARUMALITO"# Unificación de las estadísticas de los barrios "ALTAVISTA" y "CABECERA ALTAVISTA"
raw_encuesta <- raw_encuesta %>% mutate(
encuesta_calidad.barrio = case_when(
encuesta_calidad.barrio == "CABECERA ALTAVISTA" ~ "ALTAVISTA",
encuesta_calidad.barrio == "CIUDADELA NUEVO OCCIDENTE" ~ "CABECERA URBANA CORREGIMIENTO SAN CRISTÓBAL",
encuesta_calidad.barrio == "AREA DE EXPANCION SAN CRISTOBAL" ~ "ÁREA DE EXPANSIÓN SAN CRISTÓBAL",
encuesta_calidad.barrio == "CABECERA SAN CRISTÓBAL" ~ "CABECERA URBANA CORREGIMIENTO SAN CRISTÓBAL",
encuesta_calidad.barrio == "SAN CRISTOBAL" ~ "CABECERA URBANA CORREGIMIENTO SAN CRISTÓBAL",
encuesta_calidad.barrio == "PROGRESO Nº 2" ~ "PROGRESO",
TRUE ~ encuesta_calidad.barrio
)
)
# One-hot encoding para variable estudio
variables_categoricas <- raw_encuesta[, c("encuesta_calidad.p_45")] %>% mutate(
encuesta_calidad.p_45 = as.factor(encuesta_calidad.p_45)
) %>% rename(
nivel_estudio = encuesta_calidad.p_45
)
variables_categoricas_encoding <- one_hot(as.data.table(variables_categoricas))
raw_encuesta <- data.frame(raw_encuesta, variables_categoricas_encoding)Creación de la base de datos
#Agregación de las variables
db_escolaridad <- raw_encuesta %>%
mutate(encuesta_calidad.barrio = str_replace(encuesta_calidad.barrio, "ANDALUCIA", "ANDALUCÍA") %>%
str_replace("Nº 2", "NO.2") %>%
str_replace("Nº 1", "NO.1") %>%
str_replace("Nº 3", "NO.3") %>%
str_replace("AREA EXPANSION", "ÁREA DE EXPANSIÓN") %>%
str_replace("EXPANCION", "EXPANSIÓN") %>%
str_replace("AREA", "ÁREA") %>%
str_replace("BOMBONA", "BOMBONÁ") %>%
str_replace("LA ASOMADERA", "ASOMADERA") %>%
str_replace("BELALCAZAR", "BELALCÁZAR") %>%
str_replace("CALAZANS", "CALASANZ") %>%
str_replace("COLON", "COLÓN") %>%
str_replace("MIRA FLORES", "MIRAFLORES") %>%
str_replace("BARRIO FACULTAD DE MINAS", "FACULTAD DE MINAS") %>%
str_replace("CABECERA SAN ANT DE PR.", "SAN ANTONIO DE PRADO") %>%
str_replace("CARLOS E RESTREPO", "CARLOS E. RESTREPO") %>%
str_replace("URQUITA", "URQUITÁ") %>%
str_replace("LOS CERROS EL VERJEL", "LOS CERROS EL VERGEL") %>%
str_replace("CAYCEDO", "CAICEDO") %>%
str_replace("VALDES", "VALDÉS") %>%
str_replace("CERRO EL VOLADOR", "B. CERRO EL VOLADOR") %>%
str_replace("MOSCU", "MOSCÚ") %>%
str_replace("JOSELA", "JOSÉ LA") %>%
str_replace("JOSE", "JOSÉ") %>%
str_replace("EL YOLOMBO", "YOLOMBO") %>%
str_replace("PIEDRAS BLANCAS", "PIEDRAS BLANCAS - MATASANO") %>%
str_replace("BASILIA", "BRASILIA") %>%
str_replace("VILLA TINA", "VILLATINA") %>%
str_replace("LILIAM", "LILLIAM") %>%
str_replace("BOLIVAR", "BOLÍVAR") %>%
str_replace("CORREGIMIENTO PALMITAS", "PALMITAS SECTOR CENTRAL") %>%
str_replace("INES", "INÉS") %>%
str_replace("FE", "FÉ") %>%
str_replace("LUCIA", "LUCÍA") %>%
str_replace("SABIO", "SAVIO") %>%
str_replace("BERMEJAL- LOS ÁLAMOS", "BERMEJAL-LOS ÁLAMOS") %>%
str_replace("BOLÍVARIANA", "BOLIVARIANA") %>%
str_replace("EL NOGAL - LOS ALMENDROS", "EL NOGAL-LOS ALMENDROS") %>%
str_replace("JUAN XXIII - LA QUIEBRA", "JUAN XXIII LA QUIEBRA") %>%
str_replace("PROGRESO Nº 2", "EL PROGRESO") %>%
str_replace("MARIA", "MARÍA") %>%
str_replace("PLAYÓN", "PLAYON") %>%
str_replace("EL SOCORRO / LA GABRIELA", "EL SOCORRO") %>%
str_replace("FÉRRINI", "FERRINI") %>%
str_replace("LA CANDE LARIA", "LA CANDELARIA") %>%
str_replace("EL PLAYON", "PLAYÓN") %>%
str_replace("IGUANA", "IGUANÁ") %>%
str_replace("MARÍA CANO - CARAMBOLAS", "MARÍA CANO-CARAMBOLAS") %>%
str_replace("DE ABURRA", "DEL ABURRÁ") %>%
str_replace("ALTAVISTA CENTRAL", "ALTAVISTA SECTOR CENTRAL") %>%
str_replace("SECTOR CENTRAL", "CENTRO ADMINISTRATIVO") %>%
str_replace("ALTAVISTA CENTRO ADMINISTRATIVO", "ALTAVISTA SECTOR CENTRAL") %>%
str_replace("SANTA ELENA CENTRO ADMINISTRATIVO", "SANTA ELENA SECTOR CENTRAL") %>%
str_replace("PALMITAS CENTRO ADMINISTRATIVO", "PALMITAS SECTOR CENTRAL") %>%
str_replace("PROGRESO", "EL PROGRESO")
) %>%
group_by(encuesta_calidad.barrio, encuesta_calidad.comuna) %>%
summarize(n = n(), edad_promedio = mean(`encuesta_calidad.p_18`, na.rm = TRUE),
proporcion_lectoescritura = sum(`encuesta_calidad.p_35` == 1, na.rm = TRUE) / sum(!is.na(`encuesta_calidad.p_35`), na.rm = TRUE),
proporcion_escolaridad_actual = sum(`encuesta_calidad.p_36` == 1, na.rm = TRUE) / sum(!is.na(`encuesta_calidad.p_36`), na.rm = TRUE),
estudio_ultimo_ano = sum(`encuesta_calidad.p_37` == 1, na.rm = TRUE) / sum(!is.na(`encuesta_calidad.p_37`), na.rm = TRUE),
proporcion_instituciones_publicas = sum(`encuesta_calidad.p_49` == 1, na.rm = TRUE) / sum(!is.na(`encuesta_calidad.p_49`), na.rm = TRUE),
proporcion_nivel_estudio_0 = sum(`nivel_estudio_0` == 1, na.rm = TRUE) / sum(!is.na(`nivel_estudio_0`), na.rm = TRUE),
proporcion_nivel_estudio_1 = sum(`nivel_estudio_1` == 1, na.rm = TRUE) / sum(!is.na(`nivel_estudio_1`), na.rm = TRUE),
proporcion_nivel_estudio_2 = sum(`nivel_estudio_2` == 1, na.rm = TRUE) / sum(!is.na(`nivel_estudio_2`), na.rm = TRUE),
proporcion_nivel_estudio_3 = sum(`nivel_estudio_3` == 1, na.rm = TRUE) / sum(!is.na(`nivel_estudio_3`), na.rm = TRUE),
proporcion_nivel_estudio_4 = sum(`nivel_estudio_4` == 1, na.rm = TRUE) / sum(!is.na(`nivel_estudio_4`), na.rm = TRUE),
proporcion_nivel_estudio_5 = sum(`nivel_estudio_5` == 1, na.rm = TRUE) / sum(!is.na(`nivel_estudio_5`), na.rm = TRUE),
proporcion_nivel_estudio_6 = sum(`nivel_estudio_6` == 1, na.rm = TRUE) / sum(!is.na(`nivel_estudio_6`), na.rm = TRUE),
proporcion_nivel_estudio_7 = sum(`nivel_estudio_7` == 1, na.rm = TRUE) / sum(!is.na(`nivel_estudio_7`), na.rm = TRUE),
proporcion_nivel_estudio_8 = sum(`nivel_estudio_8` == 1, na.rm = TRUE) / sum(!is.na(`nivel_estudio_8`), na.rm = TRUE),
proporcion_nivel_estudio_9 = sum(`nivel_estudio_9` == 1, na.rm = TRUE) / sum(!is.na(`nivel_estudio_9`), na.rm = TRUE),
proporcion_nivel_estudio_10 = sum(`nivel_estudio_10` == 1, na.rm = TRUE) / sum(!is.na(`nivel_estudio_10`), na.rm = TRUE)
) %>% as_tibble()
# Eliminamos el barrio desconocido
db_escolaridad <- db_escolaridad %>%
filter(encuesta_calidad.barrio != "DESCONOCIDO")Escritura de los datos
write_excel_csv2(db_escolaridad, "databases/db_escolaridad.csv")Descripcion de las variables seleccionadas
- edad_promedio: se calcula como el promedio de los años cumplidos de los encuestados calculando dicho promedio para cada barrio.
- proporcion_lectoescritura: esta variable se construye a partir de la pregunta 35. La variable se calcula como la proporción de personas que pueden leer y escribir más de un parrafo con respecto al total de personas que respondieron dicha pregunta. Esta proporción se calcula a nivel de cada barrio.
- proporcion_escolaridad_actual: esta variable es un derivado de la pregunta 36. Esta variable representa la proporción de personas que están estudiando en el momento respecto al total de personas que respondieron esta pregunta para cada barrio.
- esudio_ultimo_ano: Se calcula como la proporción de personas que estudian en el ultimo año respecto al total de personas que respondieron la pregunta 37 de manera efectiva.
- ultimo_nivel_estudio: representa la moda de cada barrio extraida sobre el nivel de estudio de los hogares encuestados. Esta pregunta es extraída mediante el cálculo de la moda sobre la pregunta 45.
- tipo_institucion_educativa: se extrae mediante el calculo de la moda de cada barrio a sobre la pregunta 49. Representa el tipo de institución educativa que predominante en la cuál los encuestados estudiaron su último curso o en la que están estudiando actualmente.
Es importante considerar que la variable “edad_promedio” no necesariamente hace referencia a la dimensión de escolaridad, sin embargo, se quiere realizar un contraste entre la variable “edad_promedio” y las demás variables para extraer información y hacer validaciones empíricas.
Analisis preliminar de los datos
Medias, cuartiles y desviaciones estándar
summary(db_escolaridad)## encuesta_calidad.barrio encuesta_calidad.comuna n
## Length:308 Length:308 Min. : 4.0
## Class :character Class :character 1st Qu.: 348.8
## Mode :character Mode :character Median : 862.5
## Mean :1073.2
## 3rd Qu.:1439.0
## Max. :8987.0
## edad_promedio proporcion_lectoescritura proporcion_escolaridad_actual
## Min. :23.00 Min. :0.6667 Min. :0.1000
## 1st Qu.:32.98 1st Qu.:0.8809 1st Qu.:0.2383
## Median :36.20 Median :0.9125 Median :0.2638
## Mean :36.88 Mean :0.9079 Mean :0.2623
## 3rd Qu.:40.69 3rd Qu.:0.9421 3rd Qu.:0.2889
## Max. :58.25 Max. :1.0000 Max. :0.5000
## estudio_ultimo_ano proporcion_instituciones_publicas
## Min. :0.00000 Min. :0.2222
## 1st Qu.:0.01565 1st Qu.:0.7326
## Median :0.02206 Median :0.8948
## Mean :0.02290 Mean :0.8169
## 3rd Qu.:0.02843 3rd Qu.:0.9442
## Max. :0.11111 Max. :1.0000
## proporcion_nivel_estudio_0 proporcion_nivel_estudio_1
## Min. :0.00000 Min. :0.00000
## 1st Qu.:0.07585 1st Qu.:0.06916
## Median :0.12006 Median :0.11016
## Mean :0.12759 Mean :0.11103
## 3rd Qu.:0.16647 3rd Qu.:0.14499
## Max. :0.45455 Max. :0.25397
## proporcion_nivel_estudio_2 proporcion_nivel_estudio_3
## Min. :0.05263 Min. :0.00000
## 1st Qu.:0.19048 1st Qu.:0.06058
## Median :0.25095 Median :0.07895
## Mean :0.23695 Mean :0.07716
## 3rd Qu.:0.28760 3rd Qu.:0.09458
## Max. :0.52381 Max. :0.25000
## proporcion_nivel_estudio_4 proporcion_nivel_estudio_5
## Min. :0.09091 Min. :0.00000
## 1st Qu.:0.20988 1st Qu.:0.03508
## Median :0.24640 Median :0.04964
## Mean :0.24372 Mean :0.04925
## 3rd Qu.:0.27755 3rd Qu.:0.06209
## Max. :0.50000 Max. :0.13333
## proporcion_nivel_estudio_6 proporcion_nivel_estudio_7
## Min. :0.00000 Min. :0.00000
## 1st Qu.:0.01568 1st Qu.:0.02220
## Median :0.03062 Median :0.05376
## Mean :0.03364 Mean :0.09698
## 3rd Qu.:0.04996 3rd Qu.:0.14323
## Max. :0.11111 Max. :0.36842
## proporcion_nivel_estudio_8 proporcion_nivel_estudio_9
## Min. :0.0000000 Min. :0.0000000
## 1st Qu.:0.0004705 1st Qu.:0.0000000
## Median :0.0041754 Median :0.0009129
## Mean :0.0161616 Mean :0.0061769
## 3rd Qu.:0.0198397 3rd Qu.:0.0064606
## Max. :0.1225364 Max. :0.0727763
## proporcion_nivel_estudio_10
## Min. :0.000000
## 1st Qu.:0.000000
## Median :0.000000
## Mean :0.001348
## 3rd Qu.:0.001084
## Max. :0.024691Desviaciones estándar
db_escolaridad %>% select_if(is.numeric) %>%
apply(2, sd)## n edad_promedio
## 1.019528e+03 5.418267e+00
## proporcion_lectoescritura proporcion_escolaridad_actual
## 4.521456e-02 4.694050e-02
## estudio_ultimo_ano proporcion_instituciones_publicas
## 1.343127e-02 1.811589e-01
## proporcion_nivel_estudio_0 proporcion_nivel_estudio_1
## 6.288825e-02 5.129729e-02
## proporcion_nivel_estudio_2 proporcion_nivel_estudio_3
## 7.227418e-02 2.751268e-02
## proporcion_nivel_estudio_4 proporcion_nivel_estudio_5
## 4.962019e-02 2.082606e-02
## proporcion_nivel_estudio_6 proporcion_nivel_estudio_7
## 2.218821e-02 9.975298e-02
## proporcion_nivel_estudio_8 proporcion_nivel_estudio_9
## 2.558032e-02 1.158639e-02
## proporcion_nivel_estudio_10
## 3.221845e-03Dada la información sobre las medias arrojada para cada variable, hay que resaltar que la proporción de personas con habilidades lecto-escritoras es alta con un 0.907902, la cual es una cifra positiva pues muestra que en la población no hay analfabetismo crítico.
Por otro lado, las proporción de personas que estudiaron en el último año tambien es baja, pero no se debe concluir que esto es necesariamente negativo. Cabe la posibilidad que varias personas en un determinado barrio sean trabajadores y por tanto no dedican mas tiempo al estudio a partir de cierto momento de su vida.
>>>>>>> CONCLUIR NIVEL DE ESTUDIO AQUI <<<<<<<
Graficas
Correlaciones y análisis de distribuciones
ggpairs(data = db_escolaridad %>%
dplyr::select(-n, -encuesta_calidad.barrio, -encuesta_calidad.comuna),
columns = 1:ncol(db_escolaridad %>%
dplyr::select(-n, -encuesta_calidad.barrio, -encuesta_calidad.comuna)),
mapping = aes(alpha = 0.5))
De la gráfica anterior, se debe resaltar significativamente las correlaciones entre las variables “edad_promedio” con las variables “proporcion_lectoescritura” y “proporcion_escolaridad_actual”. Note que la “edad_promedio” esta correlacionada de manera directa con la variable “proporcion_lectoescritura” y de manera inversa con la variable “proporcion_escolaridad_actual”. Esta información que se obtiene coincide con el conocimiento epírico, pues por un lado, las personas con mayor edad tienen mas probabilidad de tener como mínimo una educación primaria, obteniendo así mas habilidades lecto-escritoras, mientras que por otro lado a medida que las personas envejecen dejan de educarse y pasan a ambientes mas laborales.
Por otro lado, note que en los barrios las personas tienen tendencia a educarse en establecimientos de índole pública debido a su fácil acceso.
Valores tipicos
ggplot(data = db_escolaridad) +
geom_histogram(mapping = aes(x = proporcion_lectoescritura), binwidth = 0.002)
ggplot(data = db_escolaridad,
mapping = aes(x = proporcion_escolaridad_actual)) +
geom_histogram(binwidth = 0.002)
ggplot(data = db_escolaridad,
mapping = aes(x = estudio_ultimo_ano)) +
geom_histogram(binwidth = 0.0005)
Escalamiento de datos
db_escolaridad_escalado <- db_escolaridad %>%
mutate_at(vars(-n, -encuesta_calidad.barrio, -encuesta_calidad.comuna), scale)Aplicación de PCA PCA AQUI <==================================
Agrupamiento mediante K-means
Selección del \(K\) óptimo (Gráfica de codo)
datos_escalados_agrupamiento <- db_escolaridad_escalado %>%
dplyr::select(-n, -encuesta_calidad.barrio, -encuesta_calidad.comuna)
k_maximo <- 15
wss <- sapply(1:k_maximo,
function(k) {
kmeans(datos_escalados_agrupamiento, k, nstart = 10, iter.max = 15 )$tot.withinss
})
plot(1:k_maximo, wss,
type="b", pch = 19, frame = FALSE,
xlab="Número de grupos [K]",
ylab="Total de suma de cuadrados intra-grupos")
k_optimo_kmeans <- 8Agrupamiento con \(K\) óptimo
agrupamiento_kmeans <- kmeans(datos_escalados_agrupamiento, centers = k_optimo_kmeans, nstart = 10, iter.max = 15 )Mapa para el agrupamiento mediante K-means
political <- shapefile("Barrio_Vereda/Barrio_Vereda.shp")
Encoding(political@data$NOMBRE) <- "UTF-8"
political$NOMBRE <- political$NOMBRE %>% toupper() %>% str_replace("DE MESA", "DE MESA")
grupos_barrios <- data.frame(barrio_nombre = db_escolaridad$encuesta_calidad.barrio, grupo = agrupamiento_kmeans$cluster)
nombres_mapa <- data.frame(nombre_barrio = political$NOMBRE)
vector_nombres = c()
vector_grupos = c()
for(nombre_mapa in nombres_mapa$nombre_barrio) {
grupo <- grupos_barrios[grupos_barrios$barrio_nombre == nombre_mapa, 2][1]
vector_nombres <- c(vector_nombres, nombre_mapa)
vector_grupos <- c(vector_grupos, grupo)
}
mapa_grupos <- data.frame(
nombre_barrio = vector_nombres,
grupo = vector_grupos
)
factpal <- colorFactor(rainbow(k_optimo_kmeans), mapa_grupos$grupo)
leaflet(data = political) %>%
addTiles() %>%
addProviderTiles(providers$OpenStreetMap) %>%
addPolygons(fill = TRUE, stroke = TRUE, weight = 2, color = ~factpal(mapa_grupos$grupo),
label = as.character(political$NOMBRE),
popup = as.character(political$NOMBRE)) %>%
addLegend("bottomright", colors = "#03F", labels = "Barrios y veredas")